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Abstract 

We study the effects of the Coulomb interaction on the growth of unstable 

modes in asymmetric nuclear matter. In order to compare with previous 

calculations we use a semiclassical approach based on the linearized Vlasov 

equation. Moreover, a quantum calculation is performed within the R.P.A.. 

The Coulomb effects are a slowing down of the growth and the occurrence 

of a minimal wave vector for the onset of the instabilities. The quantum 

corrections cause a further decrease of the growth rates. 
PACS number(s): 21.65.+f, 21.60.Ev, 21.60.Jz, 25.70.Pq 
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Very recently there has been a growing number of studies on the properties of highly 
asymmetric nuclear matter. This can be mainly ascribed to recent experimental investiga- 
tions on processes involving neutron-rich nuclei 

Under certain kinematic conditions, a nuclear collision can produce physical situations in 
which bulk or surface instabilities take place. Several theoretical models have been proposed 
in order to study such situations pj-Q. Among microscopic approaches, that based on the 
Boltzmann-Langevin equation is one of the most suitable to describe the large fluctuations 
observed in fragment formation 1^,^. In this scheme the instabilities of the self-consistent 
mean field can play a crucial role in enhancing mass and charge fluctuations 0. 

In reactions between neutron-rich nuclei regions of high charge asymmetry, like the neck 
zone in semi-peripheral collisions, can be produced. In these regions the chemical ( or 
diffusive ) instability can be the most effective 0. 

The thermodynamics of phase separation in asymmetric nuclear matter was first treated 
in studies concerning the structure of neutron stars [0. More recently thermodynamic 
approaches have been developed with specific attention to critical situations that can be 
reached in collisions of heavy ions P,pil| . 

From thermodynamics one can obtain both the equilibrium configurations and the condi- 
tions for the onset of the chemical instability, i. e. the critical values of neutron and proton 
densities and of temperature for a given equation of state. Thus it is possible to draw the 
borders of the instability in the pressure-density plane for instance. The result is that for 
asymmetric nuclear matter the region of mechanical instability is embedded in the region of 
chemical instability [§. This means that in asymmetric nuclear matter the effective spinodal 
region is defined by the chemical instability. 

In this report, rather than with the thermodynamics of phase separation, we are con- 
cerned with the dynamical development of the instabilities in asymmetric nuclear matter. 
The growth rates of the various unstable modes in asymmetric nuclear matter have been 
already studied in Refs. |n],[T2| using a Skyrme force for the effective nucleon-nucleon in- 



teraction. Here we want to evaluate also the effects of the Coulomb force on the growth of 



instabilities. It is well known that the Coulomb force gives a divergence in the energy density 
of infinite nuclear matter. Thus it cannot be taken into account in thermodynamic studies 
of such a physical system. However we will see that within a mean-field treatment of nu- 
clear excitations, the contribution from the Coulomb force does not give rise to divergences. 
Moreover, we assume that the unperturbed initial state is uniform and homogeneous. In 
this case the mean field of the initial state does not appear in calculations explicitly. Only 
the temperature, the density and the asymmetry of the initial state actually occur in the 
relevant equations. 

In the study of the nuclear dynamics by means of kinetic equations, the self-consistent 
mean field is usually treated in semiclassical approximation . In order to compare with 
previous calculations, here we study the unstable modes of asymmetric nuclear matter by 
using the linearized Vlasov equation. This equation is a semiclassical approximation of the 
RPA, valid in the long-wavelenght limit. However, we also evaluate the quantum corrections 
to the semiclassical results. 

Since the RPA is a linear approximation, we expect that the validity of these calculations 
is limited to times close enough to the onset of instabilities. An estimate of the time interval, 
in which the RPA can be considered valid, is given in Ref. |]T2||. There, the numerical solutions 
of the non linear Vlasov equation have been compared to the analytical results of the linear 
approximation. The two procedures give quite similar results for times t ^ 150/m/c. 

In our calculations we have used a Skyrme-like form for the nuclear part of the nucleon- 
nucleon effective interaction. We start from the following simplified functional for the density 
of the nuclear potential energy: 

where g = + is the total density, qi and §2 are the proton and neutron densities 
respectively, Qeq is the density of symmetric nuclear matter at saturation and Qa = Q2 — Qi- 
For the parameters in Eq. (|1|) we take the values: 

A = -356.8MeV, B = 303.9 MeV, a=-, 
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C = 32MeV, D = 130 MeV ■ fm^, D' = 40 MeV ■ f 



m 
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The values of A, B and a reproduce the binding energy ( 15.75 MeV ) and give an 
incompressibihty modulus of 201 MeV for symmetric nuclear matter at saturation with 
Qeq = 0.16 /m~^. For the values of D and D' we follow the prescriptions of Ref. [jl3| and 
Ref. respectively. Finally this value of the parameter C gives a symmetry energy co- 
efficient of 28 MeV in the Bethe-Weizsacker mass formula. The energy density of Eq. (|I]) 
coincides with that used in Ref. [|l^, except for the symmetry energy term for which we 



prefer the simplest required form. 

Concerning the Coulomb energy density, we use the expression given by the Hartree-Fock 
approximation, with the Fock term evaluated in the local density approximation: 

^-'(r)4..W/<^.'^-^©V.f. (2) 

Within the density functional theory the potential energy functional 

E^'"'\qu02] = /rfr(f(^)+^(^)) (3) 

is the quantity, which allows to define an effective nucleon-nucleon interaction for R.P.A. 
calculations |T^. The interaction is given by the double functional derivative of E^'''°*^[gi, §2] 



(4) 



calculated in the unperturbed initial state. As expected, the Hartree term of the Coulomb 
energy yields the bare Coulomb force, whereas the remaining terms of Eq. (0) and the 
exchange Coulomb term give rise to zero-range forces. 

The R.P.A. equations for the oscillations of the proton and neutron densities, Sgi{r,t) 
and 5^?2(r, t), around their unperturbed values, can be obtained by standard procedures 
( e.g. see Ref. [|T^ ). Here we report the resulting equations for their space and time 
Fourier transforms. For a wave of frequency uj and wave-vector k, the following two coupled 
equations are obtained 
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- ^^P{k,u)(Fv + k^D - D') - —)s02ik,u) = , (5a) 
- ik,uj)(Fv + P (D - D')- —)6gi{k,uj) 

V Qeq/ 

+ (^1 - ik,uj)[Fv + e{D + D') + ^)yg2ik,uj) = , (5b) 

where 

Fy= — + {a + l)B-^ (6) 

0eq Qeq 

is the volume term of the effective nuclear interaction and 

Vi^^=^^e^9^ (7) 

is the exchange term of the Coulomb interaction. The quantities and ^oi denote the total 
and proton densities in the unperturbed initial state. In Eqs. (|) (^f\k,uj) represents the 
free particle-hole propagator, which in the long-wavelenght limit takes the form 

^T\Ku:)= - [de,^p,k,co)^ , (8) 



dep 



where 



is the mean occupation number of protons or neutrons in the state with kinetic energy 
ep = p^/2m and 

$(p, k,u;)=-^J d% , (10) 

V = p/m is the nucleon velocity. We use units such that h = 1, c = 1. The effective 
chemical potential /2j in the distribution rii, is the chemical potential measured with respect 
to the uniform mean field acting on the nucleons in the unperturbed initial state. For 
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a given temperature T it is completely determined by the nucleon densities. Therefore 
the ingredients of Eqs. (^) related to the unperturbed initial state are only the neutron 
and proton densities, besides the temperature. Moreover, we remark that the product 
(^^i\k,uj)4:Txe'^ /k"^ does not diverge when 0. 

The dispersion relation uo = uj{k) is obtained by equating to zero the determinant of the 
set of equations (||). We look for solutions of Eqs. © for values of qq, T and the asymmetry 
parameter a = {go2 — Qoi)/qo inside the region of instability f^. In general we find that the 
frequency uj{k) is or purely imaginary or purely real. The growth rate of the instability is 
given by T{k) = Imuj{k). 

In Figs. and ^ the calculated growth rate is displayed as a function of k for = 0.4 g^q 
and for various values of the parameters a and T. Figures |l| and ^ show that the Coulomb 
force causes an overall decrease of the growth rate r(A;). This decrease depends only sligthly 
on the temperature( ~ 10% ) and it is almost independent of the asymmetry. Moreover we 
observe that, when the Coulomb force is included, the wave vector k must exceed a certain 
value kmin in order to get solutions of Eqs. with r(A;) ^ 0. Below k^m the solutions 
of Eqs. (|^) correspond to undamped plasmon-like oscillations, which practically involve the 
proton density alone. In the present case k^in is about 1/5 of the Fermi momentum of the 
symmetric nuclear matter ( a = ). These two effects ( decrease of r(A;) and appearance 
of kmin ) are due to the competition between the Coulomb and nuclear forces ( responsible 
for the chemical instability ). The Coulomb force pushes the protons towards regions of 
lower density, the nuclear forces instead push the neutrons in the same direction. Explicit 
calculations, with and without the Coulomb force, confirm that the proton and neutron 
densities oscillate according to this picture. 

So far we have presented results of calculations performed in the semiclassical limit 
of R.P.A.. For infinite systems this limit essentially coincides with the Vlasov equation. 
However, Figs. ^ and ^ show that the values of the wave vector of the most unstable modes 
is not a small fraction of the proton and neutron Fermi momenta ( for a = 0, k/pp ~ 0.6 ). 
Therefore, an evaluation of the quantum corrections to the Vlasov equation is of interest. 
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Moreover, since for increasing values of a the Fermi momentum of the protons decreases 
appreciably, the importance of the quantum corrections can depend on the asymmetry degree 
of the nuclear matter. For infinite systems and for the nucleon-nucleon interaction used 
here, it is a rather simple matter to perform full quantum RPA calculation. It is sufficient 
to substitute in Eqs. the semiclassical version of the free particle-hole propagator given 
by Eq. (§) with its full quantum expression. This can be obtained by a straightforward 



generalization to finite temperature of the expression for T = given in Ref. ||T^ . 

In Figs. 1^ and ^ we show the effects of the quantum corrections on the growth rates of 
the unstable modes. We observe a further quenching of the dispersion curves r(A;). To be 
more specific, the maxima of r(A;) are lowered by a factor of ~ 0.85 for all the values of 
temperature and asymmetry parameter a shown. 

For symmetric nuclear matter, a similar evaluation of quantum effects has been done 



in Ref. [|T^, without the Coulomb interaction and with a different nuclear interaction. Ne- 
glecting the Coulomb interaction, for a = we obtain results in qualitative agreement with 
those of Ref. [|T7|. 



In conclusion, for the growth rate of the unstable modes the net result of the Coulomb 
effects and of the quantum corrections can be summarized by a suppression factor of ~ 0.7, 
which does not significantly change in the range of temperature and asymmetry considered 
here. Moreover, the long range character of the Coulomb force inhibits the formation of 
chemical instabilities of wavelength A > ^-ajkmin ( — S2 fm for qq = OAQf^g ). 
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FIGURES 

FIG. 1. Growth rates of the unstable modes calculated in the semiclassical approximation 
^ Eqs. dll) ^ for T = and for three different values of the asymmetry parameter a. From top to 
bottom a. = 0.0, 0.3, 0.6. The density of nuclear matter is qq = OAQeq . Results with the complete 
interaction ( solid lines ) and without the Coulomb interaction ( dashed lines ). 

FIG. 2. The same as Fig. |, but for T = 5MeV 

FIG. 3. Growth rates of the unstable modes, calculated with the quantum corrections ( solid 
lines ) and in the semiclassical approximation ( dashed lines ). The values of the parameters are 
the same as in Fig. ||. 



FIG. 4. The same as Fig. |, but for T = 5 MeV. 
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Fig, 1 
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